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RESEARCH PROGRAM PLAN 


OBJECTIVES 

1. Technological assessment of ribbon growth of silicon 
by a capillary action shaping technique. 

2. Economic evaluation of ribbon silicon grown by a 
capillary action shaping technique as low-cost silicon. 

SYNOPSIS OF PROGRAM OF STUDY 

1. Crystal growth of silicon ribbons. 

2. Characterization of silicon ribbons. 

3. Economic evaluations and computer-aided simulation 
of ribbon growth. 



FOURTH QUARTER HIGHLIGHTS 


o A new capillary action shaping technique die was 
designed. This new die allows improvements in surface 
smoothness and in SiC surface-particle density. 


o 


Forty-seven ribbons greater than 0.5 meter long and 
25 mm wide were grown during the last quarter. 


Ribbon width was extended to 38 mm (1 1/2 inches) 


o Surface films on ribbon surfaces were analyzed as SiC 
crystallites. Significant structural differences, 

depending on the deposition location, were found. 

o Epitaxial growth of SiC through preferential 


incorporation 

of 

(111) 

SiC planes 

parallel 

to 

(111) silicon 

planes 

was 

identified as 

an important 

mechanism for 

surface 

film 

formation . 




o Development of a new technology-forecasting technique 
was continued. This technique is being applied to 
projecting the future cost of energy at the level of 
silicon-sheet material. From a baseline, future 

technology capability is projected through full 


v 



maturity. The concept of chronology is introduced by 
estimating the probability of meeting the objective 
associated with the production-unit parameter and leads 
to a specific cos t-versus- time relationship. 

o Silicon-sheet technology was shown to have the 
potential for achieving future low-cost material 
objectives for photovoltaic applications, if 
development milestones defined in this study are met. 

o 1980 and 1 985 energy-capacity costs of $750/kWE 
and $350 /kWE , respectively, at the level of 
silicon-sheet material, are projected. 

o This analysis confirmed, from a silicon-sheet material 
standpoint, that ERDA-stated energy-capacity cost 
objectives at the array level are achievable, but with 
little margin for error. 


VI 



CRYSTAL GROWTH 


by 

T. F. Ciszek 
1. INTRODUCTION 

v 

The crystal-growth method under Investigation is a capillary 
action shaping technique. Meniscus shaping for the desired 
ribbon geometry occurs at the vertex of a vettable die. As 
ribbon growth depletes the melt meniscus, capillary action 
supplies replacement material. The configuration of the 
technique used in our initial studies is shown in Fig. 1 and 
is similar to the edge-defined, film-fed growth (EFG) 

process described by LaBelle (1). The crystal-growth method 
has been applied to silicon ribbons for several 
years (2,3,4), and long ribbons up to 25 mm in width have 
been produced. 

Certain problems still await solution before the technique 
becomes viable for large-scale economical photovoltaic 
applications. High-density graphite fulfills the durability 
and wettability requirements of a die (2) and has been used, 
to date, for most silicon ribbon growth; it is not, however, 

i ■ ■ 

completely non-reac tive . Good crystallographic perfection 
has been achieved on small ribbon segments (2,3), but the 
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Crucible 


structure of large ribbons is marred by planar, line, and 
point defects. Our objective in this work is to attain a 
clear technological assessment of silicon ribbon growth by 
the capillary action shaping technique and to enhance the 
applicability of the technique to photovoltaic power device 
material . 

In this report, a new capillary die design is described. 
It represents a departure from the die types used for 
edge-defined, film-fed growth, in that the bounding edges of 
the die top are not parallel or concentric with the growing 
ribbon. The new dies allow a higher melt meniscus with 
concomitant improvements in surface smoothness and freedom 
from SiC surface particles, which can degrade perfection. 


Also in 

this reporting period. 

ribbons were grown 

for 

delivery 

to JPL. Twenty ribbons 

and 30 ribbon samples 

were 

shipped . 

Detailed dimensional 

characteristics of 

most 

r ibbons 

grown during the past 

year are presented in 

Appendix 

1 . 



Finally , 

our initial progress 

in the growth of 

38-mm 


(1-1 / 2-inch ) -wide ribbons, up to 46 cm in length, is 
reported . 
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2. CAPILLARY ACTION SHAPING TECHNIQUE 
DIE DESIGN CONSIDERATIONS 

The capillary die design used for meniscus shaping In our 
capillary action shaping technique of ribbon growth over the 
past several months Is Indicated In Fig. 2. The die is a 
departure from the edge-defined, film-fed technique of 
crystal growth. In that the cross section of the growing 
ribbon is not concentric with the top edges of the die. The 
die top is considerably thicker in the mid-region than at 
the edges, yet the resultant ribbon is flat or somewhat 


f 
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thicker at the edges than in the middle. This die design 

» 

allows a higher meniscus at the central region of the 
solidification front and thus reduces problems which can 
occur when using flat-top dies or curved-top dies with 
parallel top edges as in the EFG technique. Such problems 
come about primarily from the fact that the graphite die 
used for silicon growth slowly dissolves in the liquid. 
Carbon-saturated silicon rises up the capillary slot in the 
die and comes to the top region where ribbon solidification 
occurs. This top region is the coolest region in the growth 
system. Here, excess carbon is forced out of the saturated 
silicon solution in the form of g-SiC crystallites, which 
tend to collect at the top surface of the die. These 

crystallites distort the melt meniscus and make the ribbon 

non-uniform in its surface smoothness. Because of the 

proximity of the freezing interface to the die top, the SIC 

particles are frequently incorporated in the ribbon, where 
they generate dislocations and other defects. 

It is advantageous to keep the Interface of the freezing 
ribbon as far as possible from the die top, and this can be 
accomplished with the structure shown in Fig. 2. The die 
top surface is curved so that it is higher at the edges than 
in the middle. In this way, if the ribbon's solid-liquid 
interface is maintained approximately planar, then the 
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interface is further from the die, at least in the central 
region. The central region is most critical for generation 
of defects in the silicon body. However, if the width of 
the die top is kept constant while the die top surface is 
curved, then a higher meniscus in the central region 
necessarily implies that the top of the meniscus is thinner 
there. This .would cause the ribbon to be very non-uniform 
in thickness from one edge to the other (i.e., much thinner 
in the middle than at the ends). Thus, not only should the 
die top curve downward from the ends toward the central 
region, but it must also become wider in the central region 
than at the ends, as in Fig. 2. The meniscus, then, has a 
wider base in the central region. The wider base, combined 
with the greater meniscus height in the central region, 
results in a more uniform thickness at the .solid-liquid 
interface. In summary, there are two things that are 
important to the die design: one is the curvature of the 
top surface, and the other is the widening of the die top in 
the central region. 

Figure 3 is a cross section at the center of the die before 
the seed crystal is applied. No through capillary is shown 
because this is the area that holds the two major portions 
of the die together. Figure 3b shows the edge condition 
with the ribbon in place, and Figure 3c is 'a vertical cross 
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Fig. 3. Cross section through die and ribbon: (a) central die top before seeding, 
(bl near edge of die during growth, and (c) near central region of die top 
during growth. 


section through the central region during full-width ribbon 
growth. It can be seen that, even though the top of the die 
is narrow at the ends and relatively wide in the middle, the 
ribbon thickness is essentially uniform because the freezing 
interface is close to the die top near the ends of the die, 
but higher above the die top near the middle. The top of 
the meniscus is about as wide as the bottom of the meniscus 
at the end areas (Fig. 3b). However, the cross section of 
the meniscus near the central point of the die and ribbon 
(Fig. 3c) tapers from a wide base to a narrow top. By 
proper choice of the curvature of the die top and the taper 
angle of the sides of the die, an optimum value for this 
variation of the width of the die top with position along 
the ■ die top can -'.be.' . ^ 
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The outer edges of the top die surface, that is, those edges 
which bound the lower portion of the melt meniscus from 
which the ribbon solidifies, can be considered to be 
determined by the intersection of a vertical truncated wedge 
with enclosed angle- 0 , truncated thickness - X e , and 
width - W, with a horizontal cylinder of radius R, The 
intersection is made essentially such that the cylindrical 
surface contains the short edges, X , of the wedge top. The 
resultant die top is that of Fig. 2, where X £ is the top 

surface thickness at the ends, X the top surface thickness 

m 

in the center, 0 the enclosed angle of the tapered wedge, 
R the radius of curvature of the top surface, W the width of 
the die, and 6 the difference in height from ends to center. 
The top of the die, thus, smoothly increases in thickness 
from X to X and decreases in height, by an amount 6, as we 
go from the die edge to the die middle. 

The objective in this design is to attain a high-melt 
meniscus in the central region, since proximity of the 
freezing interface to the die top is detrimental to ribbon 
perfection and surface smoothness, while still maintaining 
the proximity at the die ends to stabilize the ribbon width. 
Furthermore, this must be achieved in a smooth transition to 
facilitate the early stages of growth from seed size to full 
width. ' 
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R and 0 are chosen to optimize the values of X -X and 6 

r me 

for successful ribbon growth. These parameters are given by 



2 tan (sin~^ W ) 
2R 


X - X - 2 6 tan 0/2. 
m e r 


For the 25-mm-wide ribbons grown recently, dies were 
constructed with R-101.6 mm, 0 *40° , and X e «0.54 mm. Thus* 
Xjjj was ^0.99 mm. The typical edge thickness and middle 
thickness of the resultant ribbons were 0.42 mm and 0.21 mm, 
respectively; the exact values were dependent upori growth 
rate. Better surface smoothness and a lower SiC 

surface-particle density (< 0.1/cm ) were seen previously 
with slightly thicker dies, where X e was 0.79 mm and R and 0 
were as above. These ribbons had a typical edge thickness 
of 0.64 mm and a middle thickness of 0.42 mm. The thickness 
of the die top in the central region was 1.37 mm in this 
case. Although these ribbons have a relatively large 

deviation from flatness (on the order of 0.1 mm), with the 
edges being thicker than the middle, they are quite smooth 

on [ a local scale compared with ribbons grown from an EFG 

. . ' .. ;• • • ■■ . '■ ' ' 

die, as was shown in Fig. 7 of the second quarterly 

j ; r V' 

report (5). Local roughness variations of 10-ym maximum 

amplitude have been achieved with the new die design, 

wheretfs roughness variations of about 50-^m maximum 
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amplitude were typical with the flat top die. Silicon 
carbide particle densities with the new die are typically 
<0.1/cm^ as compared with 5/cm^ for the EFG die. 

During this quarter, 25-mm-wide dies with the dimensions 
indicated above were used in conjunction with our standard 
growth setup [see pp.7-9 of the third quarterly report (6)] 
to grow ribbons for delivery to JPL . Forty-seven ribbons, 
corresponding to a total length of 24 m, were grown. The 
average length was 0.51 m, the average edge thickness was 
0.41 mm, and the average central thickness was 0.27 mm. 
Twenty complete ribbons and 30 ribbon segments were shipped 
to JPL. Ribbons as thin as 0.30 mm at the edge and 0.10 mm 
in the central region were produced when pull speeds in 
excess of 3 cm/min were employed. However, these ribbons 
exhibited a non-flat surface, with undulating bulges of 
small amplitude. In general, surface roughness, SiC 
particle density, asymmetrical growth, and growth-control 
requirements were all more severe with dies designed for 
thin ribbon growth than with dies designed for thicker 
ribbon growth 

Twelve^ meters of ribbon were grown from a single die in 9 
melt-down cycles. The die still appeared to be serviceable 
at this point, although ribbon roughness had increased 
somewhat because of SIC buildup at the die top. 
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3. 


GROWTH OF 38-MM-WIDE RIBBONS 

The concepts discussed relative to 25-mm-wide capillary 
action shaping technique dies were extended to the design of 
a die for 38-mm (l-l/2-inch)-wide ribbons. In this die, 
<S ■ 1.0 mm, <4 * 50° , X «* 0.44, and X " 1.16 mm. Three 
full-width ribbons were grown, the longest of which was 
46 cm at full (38-mm) width (see Fig. 4). The ribbons were 
grown at speeds of 16-18 mm/min. Typically, the edge 
thickness was 0.42 mm and the middle thickness was 0.54 mm. 
A thinner portion (0.37 mm) was present, however, between 
the edge and the middle. Thus the total deviation from 
flatness was about 0.09 mm. More difficulty was noticed 
with freeze-out during the process of spreading from 
seed-width (3-mm) to full-width growth than had been 
observed during 25-mm-wide ribbon growth. Dies are being 
fabricated with slightly different dimensional parameters, 
in hopes of reducing this problem. 

The crystallographic defect structure and surface smoothness 
of the 38-mm-wide ribbons are similar to those seen for 

25-mm-wide ribbons. The largest . ribbon grown 

(3.8 cm x 46 cm) had a SiC surface-particle density of 
0.07/cm^, which is also comparable to densities seen with 
25-mm-wide ribbons . 
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4 cm 

Fig. 4. A 38-mm-wide by 46-cm-long silicon ribbon 

grown by the capillary action shaping technique. 
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Appendix 1 

DIMENSIONAL CHARACTERISTICS OF SOME RIBBONS GROWN DURING 
AUGUST 1975 - MAY 1976 


Seed Usable Thickness (mm) 


Ribbon 

Orientation 

Length 

Width (mm) 

Seed 


Tail 


Run No. 

Axis Surface 

(CM) 

Max 

Min 

Max 

Min 

Max 

Min 

5 08 05 

0 

0 

38 

8 . 6 

6 . 7 

1 . o 1 

.00 

1.03 

.00 

5 0806 

0 

0 

65 

8 . 1 

6 . 2 

1 .00 

.00 

1 . 05 

.00 

50807 

0 

0 

25 

8 . 7 

7.6 

1 .10 

.00 

1 .15 

.00 

5 0808 

0 

0 

34 

8 . 4 

8.1 

1 .06 

.00 

1 .21 

.00 

50809 

0 

0 

114 

7 . 3 

6 . 7 

1 .01 

.00 

1 .05 

.00 

50811 

0 

0 

28 

8 . 0 

7 . 1 

1 .35 

.00 

1 .35 

.00 

50 825 

0 

0 

32 

4 . 7 

4.2 

1 .58 

.00 

1 .60 

.00 

50826 

0 

0 

32 

5 . 5 

3 . 8 

1.25 

.00 

1.43 

.00 

50909 

. 0 

0 

23 

26.2 

25.8 

.70 

.35 

.67 

.39 

5 0 911 

110 

111 

24 

25.7 

24.4 

.65 

.35 

.62 

.41 

50916 

110 

100 

46 

26.0 

. 0 

.73 

.54 

.62 

. 3 3 

50919 

110 

100 

90 

26 . 0 

.0 

.55 

.23 

.65 

.36 

50926 

110 

100 

18 

26 . 0 

. 0 

.65 

.57 

.65 

.56 

5 0927 

110 

100 

40 

26 . 0 

. 0 

.53 

. 4 0 

. 54 

.34 

50928 

110 

100 

4 5 

26.0 

..0 

. 49 

.32 

... 50 

.31 

510 0 3 

110 

112 

49 

25.6 

23 . 6 

. 67 

.51 

.63 

.46 

51005 . 

110 

112 

50 

26.1 

24.7 

.63 

.44 

.66 

.45 

51007 

110 

112 

51 

26 . 4 

23 . 0 

.75 

.53 

.73 

.45 

51008 

110 

112 

44 

25 . 6 

2 3 . 7 

.63 

.40 

.66 

.43 

'51012 

100 

110 

91 

25.9 

24.3 

.69 

.49 

.68 

.47 

5 1 0 T 3 

100 

110 

20 

26.0 

25 . 0 

.64 

.42 

.68 

.51 

51014 

100 

110 

54 

25 . 7 

.21.2 

.59 

. 6 2 

.62 

.34 

51015 

110 

100 

84 

25 .1 

23.4 

.61 

.35 

.64 

.34 

51016 

110 

100 

52 

24 . 8 

17 . 0 

.73 

.48 

. 60 

.39 

51017 

110 

100 

22 

24. 7 

2 3.8 

.64 

.45 

.44 

.26 

51018 

110 

100 

23 

25 . 1 

25 .0 

.72 

.43 

.67 

.41 

51019 

110 

1 0 0 

5 6 

25.0 

23 . 6 

.67 

.44 

. 55 

. 35 

51021 

110 

100 

■ 54 : 

26.2 

24.4 

.65 

.49 

.66 

.38 

51022 

110 

10 0 

6 3;:..;";, 

25.6 

23 . 8 

.63 :.V;, 

.38 

.59 

.28 

5102 3 

110 

100 

63 

25 . 5 

24 . 1 

.62 

.35 

.59 

.39 

5102 4 

110 

100 

24 

26.2 

12.5 

.63 

.35 

.64 

.26 

51025 

110 

100 

5 8 

2 5.7 

23. 3 

.54 

.31 

. 56 

. 22 

5 1026 : 

110 

10 0 

62 

25 . 4 

2 4.6 

. 62 

.42 

.62 

. 3 7 

51027 

110 

1 0 0 

' 59 V :.- 

24 . 9 

24 . 5 


.29 

.52 

. 2 6 

51102 

110 

112 


ik -:V. V 0 

. 0 

.00 

.00 

.00 

.00 

51104 

110 

100 

59 

25.3 

■23 . 5 


.40 

.51 

.36 
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Ribbon 

Seed 

Orientation 

Usable 

Length 

Width (mm) 

Thickness (mm) 
Seed Tail 


Run No. 

Axis 

Surface 



Max 

Min 

Max 

Min 

Max 

Min 

51107 

110 

112 

49 

25 . 5 

24.1 

,45 

. 29 

.50 

.28 

51110 

110 

100 

10 

25.7 

23 . 3 

.55 

.50 

.55 

.50 

51114 

110 

100 

48 

25.4 

24.4 

.50 

. 36 

.50 

.34 

51115 

110 

100 

87 

25 . 6 

25 . 2 

.46 

.30 

.57 

.35 

51118 

110 

100 

57 

25 . 2 

2 4.8 

.46 

. 31 

. 45 

.29 

51119 

110 

100 

54 

24 . 7 

24 . 5 

.43 

.28 

.46 

.29 

5112 0 

110 

100 

21 

24 . 9 

23.8 

.45 

.32 

. 47 

.25 

51201- 

110 

100 

16 

24 . 6 

17 . 4 

.47 

,50 

.50 

.35 

51202 

110 

100 

60 

25.0 

24.1 

. 52 

.42 

.50 

.45 

51203 

110 

100 

19 

25 . 5 

24 . 1 

.76 

.50 

.45 

.36 

512 0 5 

110 

100 

3 6 y: 

25 . 0 

24.9 

,55 

. 45 

.50 

, 45 

51207 

110 

100 

20 

24 . 5 

22 , 2 

.61 

.40 

.26 

.22 

51208 

110 

100 

46 

24, 6 

24. 3 

. 50 

. 40 

. 46 

. 3 6 

51209 

110 

112 

53 

. 0 

, 0 

.45 

.35 

. 0 0 

.00 

51210 

110 

112 

4 6 

24. 9 

20 . 6 

.67 

.50 

.65 

.50 

5 1,2 1 3 

110 

112 

32 

23 . 5 

14.0 

.38 

. 32 

,53 

.40 

51214 

110 

100 

5 5 

24 . 8 

23.2 

.48 

. 40 

. 45 

.41 

51216 

110 

100 

3 9 

25.3 

22 . 7 

. 71 

,52 

.73 

. 45 

51217 

110 

10 0 

57 

24 . 7 

23 . 4 

,74 

.48 

.75 

.47 

51218 

110 

100 

45 

22 . 9 

22.4 

.55 

.50 

. 60 

.50 

6010 2 

110 

100 

56 

25 . 5 

24. 9 

.53 

.39 

.50 

.33 

60103 

110 

10 0 

51 

25 . 2 

24 . 2 

. 48 

.31 

. 45 

. 27 

60104 

110 

112 

6 0 

25.4 

24 . 3 

.49 

. 29 

.52 

. 33 

60105 

110 

100 

58 

24 . 9 

18 . 8 

.50 

.45 

.50 

. 45 

60110 

110 

100 

46 

24 . 8 

23 . 6 

.53 

. 42 

.50 

.37 

60111 

110 

100 

94 

25.4 

23 . 3 

. 00 

. 0 0 

. 00 

. 00 

6 0112 

110 

100 

9 0 

24. 9 

23.4 

,45 

. 35 

.50 

.27 

60114 

110 

100 

33 

25.4 

22 . 7 

.50 

.45 

■ .53 ,:;t 

. 3 6 

60115 

110 

100 

36 

25 .2 

20 .1 

. 50 

. 40 

. 40 

. 3 0 

60116 

110 

100 

31 

25 . 2 

22 . 4 

.53 

. 42 

.37 

.28 

60119 

110 

100 

28 

25.4 

25 . 1 

.55 

.50 

.54 

.38 

60120 

110 

100 

9 5 

. 0 

.0 

.00 

.00 

.00 

.00 

6020 2 

110 

100 

96 

24.6 

2 3.7 

.50 

.32 

.47 

.30 

6 020 3 

110 

100 

25 

25.5 

22 . 2 

- ,57 

.38 

. 55 

.25 

6 0206 

110 

10 0 

9 

23.1 

18 . 1 

,62 

.70 

.60 

.64 

6 0209 

110 

10 0 

45 

25 .9 

24 . 2 

.80 

.90 

.64 

.90 

60212 

110 

111 

/ 5 5 

2 5 ,7 

25.1 

.60 

■ .50 

.60 

. 40 

60213 

110 

111 

77 

25 . 6 

25 . 4 

, 60 

.85 

o 

o 

.00 

60216 

110 

10 0 

40 

24. 6 

24.2 

.70 

.80 

.55 

.65 

6 0217 

110' 

100 

1 8 

24.4 

2 4 . 2 

. 60 

.85 

.42 

.52 

6 0 3 03 

110 

111 

55 . 

24. 9 

2 2 .5 

* .47 

. 43 

.45 

,35 

6 0 3 0 4 

10 0 

110 

55 

25 . 5 

22 . 6 

.45 

.42 

.42 

.22 

6 0309 

10 0 

110 

53 

25.6 

23 . 4 

.47 

.50 

Y . 42 

. 45 

6 0 31 0 

100 

.110 . 

23 

25 . 2 

23 . 9 

.35 

. 35 

.38 

.40 

60312 

111 

110 

41 

24. 7 

21.7 

.32 

.32 

... .35 y 

. 37 

60318 

10 0 

110 

4 4 

25. 6 

23.7 

.44 

. 40 

.45 

.37 

6 0320 

100 

110 

• . 5 5 

25 .3 

24 . 3 

- .42 

.38 

.36 

, 3 0 

60322 1 

111 

■112 

57 

24.7 

24 . 0 

.40 

.34 

. 41 

. 25 
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Seed 


Ribbon 

Orientation 

Run No. 

Axis 

Surface 

60323 

111 

112 

6 0 3 2 4 

111 

112 

60325 

110 

100 

603 26 

110 

100 

60328 

110 

100 

60329 

110 

112 

60330 

110 

112 

60337 

100 

110 

6 0 341 

100 

110 

60401 

100 

110 

6 04 02 

100 

110 

60404 

110 

100 

60406 

10 0 

110 

60407 

100 

110 

60408 

100 

110 

60409 

100 

110 

60414 

110 

100 

60416 

100 

110 

60418 

121 

111 

60419 

121 

111 

60420 

121 

111 

60421 

121 

101 

60422 

121 

101 

6043 2 

100 

110 

60436 

10 0 

110 

6 04 37 

121 

111 

6 043 8 

121 

111 

604 3 9 

121 

111 

60440 

121 

111 

6 0442 

121 

111 

6044 3 

121 

111 

6 0 444 

121 

111 

60446 

100 

110 

6 0.4 4,7 

121 

111 

60448 

121 

101 

60506 

100 

110 

605 11 

110 

112 

6 0 5 2 0 

110 

112 

60521 

121 

101 

605 22 

110 

112 

6 0523 

110 

112’ 

6 0525 

121 

101 

6 0526 

110 

112 

605 27 

110 

112 


Usable 

Length Width (mm) 
(cm) Max Min 

29 24.6 22.0 

54 24.4 22.2 

18 23.9 18.6 

38 25.5 21.9 

35 25.3 25.0 

56 24.8 22.7 

43 25.4 24.7 

52 25.2 18.0 

55 25.1 23.7 

56 26.2 25.8 

57 25 . 7 25.2 

76 26.4 24.6 

56 25.8 25.0 

56 .0 .0 

56 25.0 24.2 

48 24.8 23.7 

56 25.6 25.0 

58 25.4 23.6 

54 25.4 20.5 

5 6 2 5.2 24.5 

32 25.2 25.0 

50 25.2 22.4 

52 25.0 23.1 

55 24.9 24.7 

52 24.8 24.1 

55 25.2 24.1 

43 25.0 22.9 

54 25.5 24.1 

56 24.5 24.3 

16 24.4 23 . 1 

38 24.3 20.5 

39 24.4 24.3 

56 24,4 17.2 

58 24.6 19.7 

53 24.4 24.1 

46 38,2 37.8 

36 39.2 37.4 

5 8 25 . 2 25 . 0 

91 24.7 22.7 

90 25.2 24.8 

63 24.7 24.5 

55 24.7 24.6 

58 24.6 23 .3 


Thickness (mm) 
Seed Tail 


Max 

Min 

Max 

Min 

. 30 

. 27 

. 32 

.22 

. 35 

.28 

,32 

. 24 

.65 

.70 

.50 

.40 

.45 

. 38 

.54 

.46 

.45 

.42 

.40 

.28 

. 30 

.25 

.40 

. 28 

o 

CO 

. 26 

.38 

. 30 

.40 

.27 

.50 

.22 

.40 

.35 

.48 

.32 

.50 

.35 

. 50 

.37 

.45 

.32 

.45 

.30 

.40 

.21 

.50 

. 32 

.52 

. 36 

.52 

. 24 

. 00 

.00 

. 00 

. 00 

. 30 

.16 

. 30 

. 15 

.30 

.15 

. 3 5 

. 10 

.45 

. 25 

. 35 

.20 

.45 

.22 

. 38 

. 1 8 

.45 

.22 

.35 

.10 

.40 

.22 

. 32 

.10 

. 32 

.19 

.35 

.15 

.48 

.37 

.48 

.32 

.45 

.25 

.43 

.25 

.40 

. 30 

. 0 0 

.00 

.48 

.41 

.00 

. 00 

. 40 

.28 

. 45 

.32 

. 35 

.23 

.28 

.18 

.45 

.29 

.42 

.30 

.40 

.25 

. 39 

.24 

.42 

. 27 

. 42 

.25 

.40 

.32 

.36 

. 16 

.42 

.19 

.45 

.20 

.38 

.25 

. 38 

.22 

.35 

.19 

.45 

.21 

.45 

.20 

.40 

.13 

,55 

.65 

.50 

.65 

.43 

.55 

.44 

. 53 

.52 

.45 

. 51 

. 45 

.40 

.38 

. 4 0 

.28 

.40 

.35 

.45 

.28 

.42 

.30 

.39 

.30 

.45 

.40 

: .45 

.40 

.40 

.20 

,40 

.25 

.40 

.20 

.50 

. 2 5 


81 25.0 24.3 

A 0.0 ENTRY INDICATES NO MEASUREMENT NAS MADE. 
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ANALYSIS OF VAPOR-DEPOSITED SILICON CARBIDE FILMS 
ON SILICON RIBBON SURFACES 
by 

K. H. Yang and G. H. Schwuttke 
1. INTRODUCTION 

The growth of perfect single-crystal silicon ribbons through 
the capillary action shaping technique by use of carbon dies 
is complicated by the formation of SiC during ribbon growth. 
As shown in the first quarterly report (1), frequently, 
small SiC crystals form in the orifice of the die. The 
carbide growth is the result of liquid transport of 
dissolved carbon from hot regions - submersed die - where 
the equilibrium carbon concentration is relatively high to 
the cold regions - die top - where it is lower. Thus a 
carbon supersaturation occurs at the die top. This 
supersaturation at the die top is enhanced through carbon 
rejection at the growth interface and is relieved through 
carbide growth. Sometimes, crystallites get detached from 
the die and float in the meniscus at the top of the die. 
Frequently, small crystals attach to the silicon ribbon 
during growth, thus destroying the perfection of the 
crystal. . 
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This report draws attention to a second mechanism operative 
in the formation of unwanted SIC crystals. This mechanism 
is based on SIC formation on the ribbon surface via vapor 

transport deposition. Vapor-phase deposition of SIC is 

I • ■ • ■ , ■ •• : \ : 

particularly active during the seeding phase and during the 

initial growth period, leading to a more or less dense SIC 
film on the ribbon surface. Such SIC films influence 
destructively the efficiency of ribbon solar cells. To 
minimize film formation it is Important to understand the 
crystallographic nature of the film as well as its mechanism 
of growth. This report relates to these problems and 
presents a complete analysis of such films. 

2. ANALYSIS OF SURFACE FILMS ON RIBBONS 

Visual inspection of seed-ribbon crystals as grown reveals 
that the seed is covered with a dull bluish film while the 
surface of the ribbon close to the interface looks dull and 
dark. The surface dullness decreases rapidly with the 
distance from the interface, and a shiny ribbon surface is 
normally obtained after 10 cm of ribbon growth. The surface 
film formation is more pronounced for lower growth speeds 
and very strong during the seeding operation because of the 
longer residence time of the seed. In the following, such 
surface films on ribbons are analyzed through optical and 
electron transmission microscopy . 
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3. 


OPTICAL MICROSCOPY OF SURFACE FILMS 


The variation of surface-film morphology with 
distance s from the seed-ribbon interface is shown in the 
photomicrographs of Figs, la to If. Figures la and lb show 
the film structure on the seed and on the ribbon, above and 
below the interface. No particular features are resolved 
optically. The film covers the silicon surface completely. 
With increasing distance from the interface, the ribbon 
surface is covered less completely and the optical 
microscope reveals we 11- developed dendritic crystal 
structures covering the ribbon surface. The number of 
dendrites on the ribbon surface decreases rapidly with 
increasing distance s from the interface. Note the 

preferential nucleation of dendrites in grain boundaries, 
shown in Figs. Id and le. Twin boundaries do not act as 
preferential nucleation sites, shown in Fig. If. Single 
isolated dendrites may form during successive ribbon growth, 
particularly if a change in growth speed occurs. 

4. TRANSMISSION ELECTRON MICROSCOPY OF SURFACE FILMS 

Opti cal microscopy cannot identify the crystallographic 
nature of the films. Therefore, a transmission electron 
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microscopy analysis was made. For this purpose, 
3-mm-diameter specimens were cut ultrasonically out of seed 
and ribbon crystals at different locations. The specimens 
were jet-etched, using a mixture of HNO 3 /HF. Before 
etching, the specimens were thinned by mechanical lapping 
(on one side only), mainly for film removal. Subsequently, 
the jet etch was applied to the lapped side until the 
specimen was thin enough for electron beam penetration. It 
was noted that the surface film was very resistant to the 
etch. Consequently, it was easy to etch holes into the 
silicon and thus isolate the film for transmission electron 
microscopy work. Figure 2 shows some typical results for 
the films covering the seed section. Figure 2a is a 
transmission electron micrograph of a seed specimen. Note 
the square structure of the silicon holes typical for the 
(ICO) orientation of the seed face. Figure 2b is the 
corresponding electron diffraction pattern, and Fig. 2c 
gives the aperture-limited dark-field pirture obtained at 
position A indicated in Fig. 2b. The diffraction pattern 
(Fig. 2b) is analyzed as 3-S1C. The results are summarized 
in Table I. The dark-field analysis (Fig. 2c) indicates 
that the film consists of randomly oriented crystallites. 
The size of these crystallites ranges approximately from 700 
to 1500' ft. ■ 
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Fig. 2. TEM micrographs showing (a) £ -SiC film on (001) seed surface, (b) electron diffraction 
pattern of (a), and (c) dark -field image taken from A in (b). 


•• * 

r- * ; 




- 




WESm. 


22 




TABLE I. Identification of Surface 

Film on Crystal Seed 

Surface Film 

d , R 

Ref lection 
HKL 

ASTM #1-119, S-SiC 

d, R 

2 .510 

111 

2 . 51 

2 . 173 

200 

2 .17 

1.541 

220 

1 . 54 

.1.2 98 

311 

1.31 

1.258 

222 

1.26 

1.089 

400 

1 . 09 

0.999 

331 

1 . 00 

0.972 

420 

0.97 

0.888 

422 

0.89 

0.8 37 

333, 511 

o 

00 

■C- 

0.767 

440 

0.77 

0.733 

531 

0.74 

Similar transmission electron micrographs of specimens taken 

from the ribbon 

at positions close to 

the interface indicate 

that the surface film on the ribbon 

also consists of SiC 

crystals . 

However, the crystallite size is 

approximately 1 

pm. Examples; of such 

crystallites are given 

in Figs. 3a- c. 

Figure 3a is a 

standard bright- field 


micrograph, while Fig. 3b is the corresponding dark-field 

' ■' • .■ ■■ "I '• ■ • ' . : . • 1 :• • ' 

picture obtained by placing the "limited" aperture at a 

section of the (111) SiC ring (Fig. 3c). This result 
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Fig. 3. TEM micrographs showing (a) the bright-field image, (b) the dark-field image of 

0 -SiC particles on ribbon surface close to seed-ribbon interface, (c) corresponding 
electron diffraction pattern. 
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identifies the grains marked A, B,- and C in Fig. 3a as 
g-SiC. Moving the aperture around the (111) SiC ring to 
different positions causes crystallites of different 
orientation to show up successively in the corresponding 
micrographs. From such observations it follows that the 
film consists of randomly oriented B-SiC crystals. 

Interesting and instructive results are obtained through the 
t r ansmi s s ion electron microscopy analy sis of the dendritic 
structures. Dendritic structures on ribbon surfaces have 
been found for growth rates of 12 mm/min to 30 mm/min. The 
dendrites occur randomly over the ribbon surface or 
preferentially along grain boundaries, and are of submicron 
size for faster growth rates. Consequently, they may not 

yield to optical inspection. The dendrites have been found 

> 

to influence generation lifetime of the silicon ribbon 
surfaces. A detailed investigation of their influence on 
generation lifetime in silicon ribbons is in progress. 

The transmission electron micrographs of a dendritic cluster 
in bright and dark fields are shown in Figs. 4a and b . The 
corresponding transmission electron diffraction pattern is 

. • ’ ■ : • ■" ' ■ = i: * ...... ■' ■ . . : . . ' . , . 

■■■■• :■■■'■ • - • j . ' • ; •' • 

presented in Fig. 5a; the diffraction pattern is reproduced 
schematically in Fig. 5b. The diffraction pattern of 
Fig. 5a contains the basic [001] silicon diffraction 
pattern, but, in addition, every silicon reflex is 
surrounded by characteristic satellite reflexes . 
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Fig. 5. (a) Selected-area electron diffraction pattern of Fig. 4, (b) schematic drawing of (a), (c) superposition 
of 0 -SiC and Si (001) diffraction patterns with Si (220) reflex p acting as a secondary source for 
double diffraction. Large and small solid dots represent Si and 0 -SiC reflections, respectively. 

The small solid dots represent the double-diffraction spots due to p. 
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For analysis of the result contained in Fig. 5a the 

diffraction pattern of Fig. 5b was reconstructed with the 
following assumptions: 

1. The dendrites consist of B-SiC and grow epitaxially 
with [001] orientation on the [001] oriented silicon. 
The epitaxial relationship between SIC and silicon is 
perfect for the (110) B-SiC planes parallel to the 
(110) silicon planes. Under such conditions, extra 
diffraction spots due to SiC appear in the silicon 
pattern. The distance between the silicon main reflex 
and any (h,k,l) diffraction spot due to B-SiC is 

obtained as follows. The ratio of lattice constants 

between B-SiC (a * 4.358 X) and silicon (a - 5.403 X) 

is 0.803. In reciprocal space, this corresponds to a 

distance of 1:0.803 ■ 1.245. Consequently, any (h,k,I) 
SiC reflection is located at the distance 1.245 <h,k,l> 
from the main beam. 


2. Additional satellite reflexes around the silicon dots 
are produced by the epitaxial B-SiC phase due to double 
reflection. Double diffraction occurs if a diffracted 
beam from the silicon passes into the epitaxial SiC, or 
vice versa. In both cases the double-diffracted beam 
is determined by adding together the reciprocal lattice 
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vectors corresponding to the two component 
diffractions. Extra reciprocal lattice points 
displaced from the silicon matrix result. 

Taking these two mechanisms into consideration and the 
additional fact that double diffraction from SiC (400) 
planes is found to be very weak, the diffraction pattern of 
Fig. 5c was constructed. The large solid spots in this 
pattern are due to silicon of <001> orientation. The small 
solid dots represent the SiC reflections, and the open-small 
dots represent the double-diffraction spots for the silicon 
reflex (220), p, acting as a secondary source. If this 
pattern is reproduced for the four (220) and four (400) 
silicon reflexes acting as secondary sources and the results 
are superimposed into a single schematic pattern, the 
superposition yields the pattern given in Fig. 5b, which 
describes exactly the diffraction pattern of Fig. 5a, 

Additional information on the epitaxial relationship between 
g-SiC and Si is obtained through tilting experiments in the 
electron microscope. Analysis of electron diffraction 
patterns of [114], [112], and [114] orientation yields the 
orientation relationships summarized in Table II . 
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TABLE II. 


Epitaxial Relationship Between $ - S i C and Si 



Orientation Relationship 

Orientation 

3-Sic 

parallel 

Si 

114 

[ 3 ll ] 


[311] 


[131] 


[131] 

112 

[ 3 1 1 ] 


[311] 


[HI] 


[111] 

114 

[311] 


[311] 


[131] 


[131] 


Further insight into the SIC growth on silicon is obtained 
from the following results. Figure 6a shows the 
bright-field micrograph of a dendrite onj a <112> ribbon 
surface. This particular surface plane is perpendicular to 
(111) crystal planes. This surface orientation is the 
result of twinning as described in the third quarterly 
report (2) . The dendrite grows again preferentially in 
<1T0> directions. This is similar to dendrites growing on 
<001> ribbon surfaces. Interesting is the dark-field image 
of this dendrite recorded through use of the { 220 ) B-SiC 


reflection, as indicated in Fig. 6c. The dark-field image 
i 8 given in Fig. 6b and shows a group of parallel layers 
spaced at approximately 500 % in the [110] direction. This 
result indicates that the dendrites consist of a succession 
of silicon and SiC layers which are stacked along (111) 



Fig. 6. TEM micrographs showing (a) bright-field image, (b) dark-field image of a dendrite on a 
< 1 12> ribbon surface, and (c) electron diffraction of (a) and (b). 
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The electron diffraction pattern of Fig. 6c indicates that 
the epitaxial relationship between g-SiC and silicon is 
(112) g-SiC parallel to (112)Si, with (220) B-SiC parallel 
to (220), and (111) 3-SiC parallel to (lll)Si. Fig. 6c is 
essentially identical with the diffraction pattern taken by 
tilting a [ 001 ] -or iented specimen into the [112] 
orientation . 


The simple epitaxial relationship observed in this study is 
in good agreement with the results reported 
previously (3-5). Brown and Watts (3) and Jacobson (4) 
reported that the growth of g-SiC on (001) silicon 
substrates, by use of chemical vapor deposition, results in 
the orientation relationship of (001) g-SiC parallel to 
( 001 ) Si with (220) g-SiC parallel to (220)Si, and (220) 
g-SiC parallel to (2l0) Si. On (110) and (111) silicon 
substrates, similar epitaxial relationships are established 
between g-SiC and the silicon substrate (3,5). 

5 SUMMARY 

Surface films on silicon ribbons! grown by the capillary 
action shaping technique by use of carbon dies are analyzed 
through optical and transmission electron microscopy. The 
films are formed through vapor deposition and consist of 
g-SiC. The SiC shows significant structural differences, 
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and on 


depending on deposition location - seed or ribbon 
ribbon growth speed. On the seed surface the SiC deposits 
as randomly oriented crystallites ranging in size from 700 
to 1500 %. Close to the seed-ribbon interface the 
crystallite size increases to 1 ym. The small crystals are 
of well-defined crystallographic shape. With the increasing 
speed of silicon ribbon growth, epitaxial formation of SIC 
dendrites on the silicon ribbon surface becomes the dominant 
SiC growth mechanism. The epitaxial growth of 8-SiC occurs 
through preferential incorporation of (lll)SiC planes 
parallel to (111) silicon planes according to the epitaxial 
relationship (001) 8-SiC parallel to (001)Si with (110) 
g-SiC parallel to (110)Si, and (llo) 0-SiC parallel to 
(lTO)Si. 

i* * 
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CURRENT OUTLOOK FOR LARGE-AREA SILICON SHEET - 
A TECHNOLOGY PROJECTION 
by 

A. Kran 

1. INTRODUCTION 

Interactive computer simulation is used to support the 
development of technological and economic data required to 
define the potential of silicon-sheet growth for large-scale 
photovoltaic applications. The silicon-ribbon growth 

production-unit model simulates the complex interactions 
between physical variables pertaining to ribbon processing 
and the economic parameters associated with product 
manufacturing and business management. 

As described in the second quarterly report (1) , the 
production-unit concept, together with technology 
forecasting and sensitivity analysis, was used to compare 
single- and multiple-ribbon growth systems for their ability 
to provide low-cost silicon-sheet material. Conclusions 
favored single-ribbon growth and suggested that 
processing-technology improvements offered the best 
potential for achieving low-cost silicon-sheet material 
objectives within the shortest period of time. This report 
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extends that work to include an assessment of 
energy-capacity cost at the sheet material level, using 
system simulation and probability concepts (2). It also 
associates technology parameters projected in time with an 
estimated probability of the event’s occurring at a specific 
point. This results in a cost ($/kWE) versus time 
relationship and in an ability to compare our projections of 
future sheet material cost with ERDA projections for 
solar-array costs. 

Since single-ribbon growth, oh the basis of previous work, 
appeared to have more potential for achieving low-cost 
material objectives, no further analysis of multi-ribbon 
growth is, planned without additional data. 

2. TECHNOLOGY FORECASTING METHODOLOGY 

2 . 1 General 

Technology forecasting, today, is recognized as an integral 
part of the decision-making process, leading to a commitment 
of resources to future products. Properly structured and 
applied, it is a useful tool for looking ahead to 
increasingly complex technology and to an environment marked 
by rapid changes. Technology forecasting is also a 
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structure for communication, requiring the forecaster to 
define his terms and open to challenge his technical 
expertise, data, thought process, and biases. It compels 

. ■; ■ ; ! ' . • I j\ ' '• ■ :/ .j 

him to make his assumptions explicit, which does not 
guarantee the correctness of forecast, but does offer the 
opportunity to confirm assumptions and to ensure that the 
data used are the best available. 


Many forecasting techniques are reported in the literature. 
Most concern themselves with the evolution of a product, 
such as a computing system, calculator, or automobile, and 
involve the fitting of some mathematical function or curve 
to historical data. The forecast is obtained by projecting 
the fitted curve into the future. Learning-curve analysis 
is an example of this type of general, aggregate projection, 


which is useful for reviewing the general viability of a 
product or technology. Once a project is under way, 
however, more specific monitoring and projection techniques 
must be added, so that progress from a "bottoms-up" point of 
view can be compared with management-stated "tops-down" 
program ob j ectives . 
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2.2 Forecasting Technique 


Our procedure for forecasting silicon-sheet manufacturing 
technology is predicated upon the production-unit concept: 

1. Define a specific production unit. 

2. Project technology capability of parameters. 

3. Introduce chronology. 

4. Evaluate numerically to obtain cost vs time. 

The production-unit approach reduces the complexity of 

t ... ■ 

interaction among processing sectors in a manufacturing 
operation, and may be thought of in terms of three elements: 
processing technology, resources, and raw materials. Its 
purpose is to transform poly crys talline bulk silicon into 
sheets of single-crystal or controlled-crystallography 
material; suitable for solar-cell fabrication. These three 
elements are described to the system as a specific 
combination of manpower, crystal-growing equipment, and 
polycrystalline silicon needed to progress through the 
crystal-pulling sector in a solar-cell manufacturing 
operation.. : 
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The "baseline" production unit, or starting point for the 
projection, is developed from specific data which reflect 
processing technology practiced in the laboratory and 
assumed transferable to a production environment, and from 
estimated direct- and indirect-cost items, plus profit, 
representative of a small*- to medium-size concern . From 
this baseline, any one of the 27 model parameters can be 
projected in terms of what is anticipated in the 
"near-future" (soon) , what is expected in the "future," and 
what will be approached as the "limit." This is done 
without commitment to a specific time frame. Subsequently, 
in order to introduce chronology, these technology 
projections are associated with a probability of meeting 
particular technology objectives by a specific point in 
time. For example (as further discussed in section 2.4), we 
estimate for 1985, with a probability of 70%, that ribbons 
between 10 and 25 cm wide will be routinely pulled. After 
numerical evaluation of the cumulative probability 
distribution function, calculation of energy cost at the 
level of silicon-sheet material, and statistical analysis of 
these cost figures, a cost-versus-time curve is plotted. 

This forecasting technique requires two projections: the 
first, relating to technology and its limits, should be 
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provided by the technologist; the second, pertaining to the 
probability of implementing stated technology objectives, is 
more management-oriented, as implementation schedules depend 
upon allocation of resources to develop the technology. 
Finally, because this technique is readily applied and 
rapidly iterated, it lends itself well as a tool to be used 
jointly by technology developers and planners interested in 
tracking development progress and cost profiles. 

2.3 Baseline Definition 

The baseline production unit may be thought of as a 
reference point from which all projections are made. The 
list of 27 input parameters to the model is shown in Fig. 1 
and is divided into three categories. The ribbon data 
category contains processing-related parameters, such as 
ribbon width (2.5 cm), growth rate (1.5 m/hr), and 
thickness (0.3 mm). Wherever possible, such as in the case 
of ribbon width, growth rate, or thickness, state-of-the-art 
values are used. 


The direct-cost category comprises the cry s tal- growth system 
cost - which, at $50,000, is essentially the cost of a 
Czochralski puller modified for ribbon growth - equipment 
life (7 years) , and interest rate •(10%) , so that equipment 
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f fflflfoN CffOUTH RATE.tl/H* 9 - 2 '?.5*>' 5 1 .50 1.50 ..05 55 PC7 

? ?{!fg N of MI Sf S M r r- feJMlW 7. 75 50 


30 30 30 


WM* , DOLLARS 


50000 50000 

7.0 7.0 7 

- 10.0 10.0 

PFffCFNf - 


50000 

0 

10.0 

70 75 


7 EQUIPMENT LIFE. YEARS 

8 INTEREST RATE. PERCENT 

9 EQUIPMENT AVAILABILITY. PERCENT - 70 75 80 

PERSONNEL PER SHIFT PER MACHINE _ 

10 11 NO. OF SUPUS ' * 0.05 0.05 0.0? AT 8 - 25000 25000 25000 

12 13 NO. OF ENGRS * 0.10 0.10 0.10 AT 8 - 20000 20000 20000 

14 15 NO. OF TECHN - 0.50 0.50 0.50 AT 8 - 10000 10000 10000 

16 POLY SILICON COST. DOLS/KG - 65 65 65 

17 POLY YIELD TO RIBBON .PERCENT - 80 83 85 

SERVICES AND SUPPLIES 

18 CRUCIBLE/DIE/PARTS COST PER DAY - 30 30 30 DOLLARS 

19 POUER COST AT - 0.05 0.05 0.05 DOLLARS PER KUH 

20 ENERGY TO OPERATE EQUIPMENT - 12 12 12 KU 


21 22 a-'H ■ 

23 - 


50 50 50 PCT OF PERS* 10 10 10 PCT OF RAN MATL COST 

10 10 10 PERCENT OF DIRECT COST ’OVERHEAD 


24 e&QELI BEFORE TAX. PERCENT - 10 10 10 OF DC’O'H’G* 


25 UORKUEEk. HOURS - ’ 168 168 168 

26 CONVERSION EFFICIENCY. PERCENT - 8.00 

27 ENERGY DENSITY AT AMI. KU/SQ M PEAK ■ 


8.00 
1 1 


6.00 

1 


Fig. 1. Baseline input parameters to production unit. 


capital recovery can be calculated. Equipment availability 
is defined as the percent of time the system is available 
for crystal pulling, excluding setup, polysilicon melt-down, 
and random machine failure time. Also included here are the 
direct personnel required to assure efficient operation of 
the system, polysilicon cost, and services and supplies, 
which include die cost. Overhead, general and 
administrative expenses, and profit are defined as 
percentages relating to other direct-cost items. 


The third category defines the workweek in terms of hours 
and of energy-conversion efficiency at AMI, a hypothetical 
value to assess energy-capacity cost at the level of 
silicon-sheet material . 


40 


Since ribbon growth is still considered to be in the 

development, not manufacturing, phase, parameters such as 
yield of "cell quality" (suitable for solar-cell 
fabrication) ribbon, poly yield to ribbon, and machine 
availability are estimated in terms of what could be 

expected in a production environment on the basis of 

laboratory experience. Since these values are subject to 

interpretation, three cases are shown, reflecting an outlook 
ranging from conservative to optimistic. 

Output from the model, shown in Table I, consists of the 

major factors contributing to sheet material and 

energy- capaci ty cost. For our purposes, the most 

conservative values (energy-capacity cost: $8476/kWE peak) 

are used as the baseline. They include the average yielded 

2 

growth rate (0,02 m /hr) , the yield factor (0.56), plus the 

2 

following direct-cost elements, calculated in dollars/m : 

equipment capital recovery, personnel, polysilicon, and 

services and supplies. Also calculated in dollars/m are 

overhead cost, G&A expense, and profit. The addition of 

2 

these items results in a total dollars/m figure for 
silicon-sheet material ($678), representing a selling price 
to a manufacturer, or purchase cost to a buyer. 
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TABLE I. Economics of Silicon Ribbon (Baseline) — 
One Ribbon Puller 


SlUUL&ILQU cr dttC ll.'04/75 

37 

38 

39 

RIBBONS GRCUN SIMUL TANEOUSL Y 

1.00 

1.00 

1.00 

RIBBON UIDTH, CM 

2.50 

2.50 

2.50 

AUG YIELDED GROUTH RTE.SQ M/HR 

0.02 

0.02 

0.02 

COMBINED YIELD FACTOR 

0.56 

0.62 

0.60 

DIRECT COST IN DOLS/SQ METER 
EQUIPMENT CAPITAL RECOVERY 
PERSONNEL 
POLY SILICON COST 
SERVICES -'SUPPLIES 

66.54 

224.49 

81.13 

71.45 

57.96 

195.56 

72.99 

63.66 

50.94 

171.08 

66.62 

57.20 

SUBTOTAL : 

443.61 

390.17 

346.64 

OVERHEAD COST IN DGLS-SQ METER 
G*A EXPENSES IN DOLS-'SQ METER 
PROFIT IN DOLLARS -'SO METER 

116.79 
56.04 
61 . 64 

102.32 

49.25 

54.17 

90.45 

43.73 

40.11 

TOTAL COST IN DOLS-'SQ METER 

678.08 

595.91 

529.16 

DOLLARS PER KU 

8476.03 

7440.89 

6614.44 


2.4 Parameter Projection and Numerical Evaluation 

All technology projections are made from the baseline and 
address three future points in time: near-future, future, 
and limit. These points define the expected capabilities of 
the particular technology, silicon-sheet growth in this 
case, from its state-of-the-art to full maturity. 
Chronology is introduced by estimating the probability of 
meeting the objective associated with the technology 
parameter at a stated point in time,. 

Ribbon-width capability is used to illustrate this 
technique. As can be seer, from Table IT, we project for 
1980, with a probability of only 0.05, that ribbons between 

2.5 cm wide, our current baseline, and 5 cm wide will be 
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TABLE H. Ribbon-Width Capability 



Present 



Near-Future 



Future 



Limit 

Width (cm) 

2.5 

4 


5 

4 


10 

y 


25 

Probability 











1980 $ 


*05 


.65 


♦ 1 

30 

■1.0 

1985 $ 


.02 


.28 


• 

♦ 

70 

= 1.0 


routinely pulled. In other words, we are convinced that 
this capability will be achieved. We estimate the same 
probability to be 0.65 that ribbon width will be between 5 
and 10 cm, and 0.30 that ribbon width will be between 10 and 
25 cm. The probabilities are mutually exclusive and add up 
to 1* For 1985, the probabilities are 0.02, 0.28, and 0.70, 
respectively. A graphical representation of- the foregoing 
is shown in Figs. 2 and 3. 


The probability is continuous and even for each point 
between the three intervals, as is seen from the 1980 
cumulative distribution function shown in Fig. 4, where we 
transform a uniformly distributed sequence of random numbers 
(y-axis) into one that is non-unif ormly distributed and that 
is based upon our estimated probability distribution 
(x-axis). This cumulative probability distribution function 
is then evaluated numerically, as seen in Figs. 5 and 6. 
For example, using the same random number (0.55) , 
Fig. 5 (1980) yields a 9-cm-wide ribbon, whereas 
Fig. 6 (1985) yields a 15 . 5-cm-wide ribbon. 

43 Technology Projection 









Fig. 5. Ribbon-width capability vs 1980 cumulative probability- 
example of numerical evaluation technique. 



Fig. 6. Ribbon-width capability vs 1985 cumulative probability- 
example of numerical evaluation technique. 


3. COMPUTER PROGRAM FOR TECHNOLOGY FORECASTING 


3.1 Data Flow Model 


The computer program is written in APL, a high-level 
programming language, and provides an interactive system for 
technology projection. Its operation is briefly described 
in conjunction with the data flow model shown in Fig. 7. 
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Fetch 
stored PU 
parameters 


Display 
list of PIT 
parameters 


i 


Calculate 



I 


y Display 
/ calculated 
\. results 



PU = production unit 


Fig. 7. Technology-forecasting data flow model. 



The program operates on a 27-by-8 matrix, containing 4 
columns of production-unit values and 3 columns of estimated 
probabilities pertaining to them (see section 4.2). As 
depicted in Fig. 7, the program fetches the first row of the 
matrix, generates a random number, develops the appropriate 
cumulative distribution function from the given 
probabilities, and then numerically evaluates the cumulative 
distribution function, using a random number generator, to 
determine the specific production-unit parameter. It then 
proceeds to the next row until a complete set (27) has been 
developed. Subsequently, energy-capacity costs (three) at 
the' silicon-sheet material level are calculated (Fig. 1) and 
stored, as are three sets of input parameters (those 
resulting in the highest and two ‘lowest $/kWE figures). 
This operation is iterated the number of times specified 
(600 cost figures each for 1980 and 1985). Finally, a 
statistical analysis provides the data points for the curve 
(shown in section 4.5). 

3.2 Listing of Functions 

Following is a display of functions needed to operate the 
program in conjunction with the code published in Ref. 1. 
The computer code is documented to facilitate understanding 
and maintenance by o thers . 
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PROJECTION OF ENERGY CAPACITY COST AT THE LEVEL OF SILICON 
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DEFINE . MODIFY . AND PROJECT . DEFINE SEIS UP T&E 
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2500 PROJECTIONS CAN BE CALCULATED , BEFORE THE UORK SPACE IS 
FULL. 


m 
(21 
[ 3 1 
1 4 3 
15 ] 
(63 

m 

161 
193 
1101 
111 ] 
{121 

( 13] 

(14] 
(15 3 
[161 
117] 
C16] 
[ 191 
[201 
[21 ) 
(221 
[231 
C241 
(251 
[26] 
[271 
[26] 
(291 
[30 1 


VDEFINEl 0 17 
7 DEFINEtAtBtCtK 

M COMPILES GLOBAL MATRIX M. 27 R0HS*4PRQJ *ZEST 
M ALSO A SAMPLE GU MATRIX FOR PU MODEL 
' USE EITHER FORMAT FOR ENTRY t< 

' 2 5 10 25 .3 .4 .3 (7 ENTRIES >’ 

1 5 10 25 ( THEN VECTOR IS XU ! 6 7H 1 1 1' 

• 2 5 10 25 3 4 3 (XC5 6 73 IS DIU BY 10)' 

A - 27 3 f0 
B+ 27 7 fO 
K+0 
LOOP i 

+ U7<K***l)lEND 
'ENTER ROU 1 ,TK 
■»<7-n>»O)t60 
oq.c,3m 
GO i 

•♦((♦/as s ?])no)tGi 

C(5 6 73KU5 6 73*10 

GU 

« ROUS 1 7 MEMBER UF.CT , */“ZrC <•> 1 OR 3 

ft 

-*< (7«ftC)*((l»*/“3tC)*(3»*/ _ 3tC) ) )tG2 
'ROU *,070,' INPUT NOT ACCEPTABLE >,T C 
K+K-l 
* LOOP 
G it 

AIK} J*PPg,a/0 ]*C 

m 

■H.OOP 
END i 


VMODIFYl 03V 
7 MODIFY XlAlK 

[1] ftft MODIFIES GLOBAL MATRIX 03 
[23 ft RIGHT ARGUMENT'" 'ROU NUMBER (S 1 
[33 ‘USE FORMAT t 2 5 7 10 .5 .3 ,2' 

[43 X*-,X 

t$] K + 0 
[6 3 LOOPi 

(71 -CO,X><*«*» INTEND 

(61 'ENTER ROU > ,r<X[KV 

[ 9 ) A-.O 

[103 ft INPUT MUST BE 7 NUMBS.' 3»«<*>1 OR 3 
till ft 

[ 12 ] ■♦( (7«ft*W(W/ 3»A).n/C0 

(131 'ROU *,(TX). ' INPUT NOT ACCEPTABLE} ' ,TA 
[14] (OX-1 

1 15 ] ♦ LOOP 

(16 3 <30 5 

(173 eStXlKl) JM 
(16 3 +LOOP 
[193 END: 

v 


48 





wwjfcnoi* 

» PROJECT XlAlBlDtK 

M CALCS PROJCTD ENERGY CAPACITY COSTS (DOLS/XUE) 
aa FROM ft LIST PROVIDED BY RNL. INPUT TO ST AT US 
a RIGHT ARGUMENT * • >NO OF SETS (3 ) TO BE EVALUATED 
a RESULT IS 3»X EN (COSTS) PLUS DISPLAY OF 

• THREE (HIGHEST, 2LOUEST) PRODUCTION UNITS 

( M . *-1000000 
EN-iB*K*0 
D - 1001000 
A + 27 3 *0 
Git 

*(X<K+K*IMPRT 
CALC 

EN+rEN,UC 
■*<K«*1)tH0 


(11 

12 ) 

(3) 

(4) 

C 5 1 
( 6 ) 
17) 
16) 
19) 
(10) 
( 11 ) 
(12) 

(13) 

(14) 

(15) 

(16) 
(17) 
( 10 ) 

(19) 

( 20 ) 
( 21 ) 
( 22 ) 

(23) 

(24) 

(25) 

(26) 
(27) 
( 20 ) 
(29) 
(31) 

(31) 

(32) 

(33) 

(34) 

(35) 

(36) 

(37) 
(36) 
(391 
(401 

(41) 

(42) 

(43) 

(44) 

(45) 
(461 
(471 


PICK MID -UAL U ft COLUMN TO INITIAL (371 


a 
a 
• 

4(;3)»cuniu*4gc] 

HO l 

• 

• !<*■> HIGHEST EN VALUE 

+(BiU/ttC))*G2 

a 

• RETAIN COLUMN OF ft RESULTING IN HIGHEST EN 

4(>l)*ftt>tfCir/aC) 

62l 

'*(Di(l/(1C))tGl 

a 

■ RETAIN ft COLMS RES IN 2 LOUEST EN VALUES 

A( J2)«*(J3] 

A(;3)«ftt itfCU/OC) 

>i/ac 
•*61 
PRTt 


sr 


INVERT 
LIST 

0 

CALC 
DISPLAY 1 


9PP*l 01» 

f Z*PPR XlAtBlK 
•• CALCULATES 


aa 

M 

* 

* 

4 

• 

a 

a 


( 1 ) 

( 2 ) 

(31 

(4) 

(5) 

(6) 

(7) 

( 0 ) 

(9) 

1 10 ) 

( 11 ) 

( 12 ) 

(13) Wit 

(14) a 

( 15) 

(161 
(17) 

(10) 

(19) 


3 PRODUCTION UNIT PARAMETERS 


FROM TECHN PROJECTION AND EST PROBABILITY, 
USING RANDOM NUMBERS 

RIGHT ARGUMENT <•> 4 TECH PROJ , 3 EST PROB 
EXPLICIT RESULT «•* 3 PROD UNIT PAR VALUES 

IF */~3tX IS 3 AND X[ 1 ) IS NOT 0. 

TECHN LIMIT HAS BEEN REACHED .PROJ IS N/A 


-*((3*“/ 3*X) V ( 0*ltX) )/Hl 
Z*-3'X( 1) 

» 


IF */‘3tX IS 3 # ' PARAMETER IS NOT BEING 
EVAL. X(l) SHOULD BE SET TO 0. 


OF POOR. QUALITY 


■*(' »/ Z*X)*Z)/H2 
Z-X12 3 4) 

-*0 

(20) H2t 

(21.1 •*((>/X(5 6 7))-l)/W3 

122 1 -*O^0*- ILITIIS MUST*!' 

(231 H3t 

(24) Z+ZH 

(25) a 

(26) a CALC SLOPE OF 3 LINE SEGMENTS .THIS IS CUM 

(27) a PROB DIST FN,Xl 11 AT V*0.X(4) AT V*1 
(201 a 

(291 »-<X(2 3 4 1-X( 1 2 3im<*\X(5 6 7))-(«N(0.X( 
5 6 )))) 


49 Technology Projection 



t30) K-0 

[ 31 1 LOOP » 

[321 • f UAL CUn PROS DIST FN (3 VALUES) 

1 33 3 • 

[34] -(3'KHTM )/0 
[ 35 ] A -Pad 
(361 • 

137] a IS RANDOM NO IN 1ST LINE SEC RANGE 7 
(301 • 

(391 '•(4’XISD/Cl 
[ 40 3 • 

[41] • PAR UAL IS X(lHim>«] 

[42] a 

[431 ZIK ]*X[ 1 ]♦(*[ 1 ]>«) 

[44] -LOOP 
(4?j an 

[46] a 

[47 1 a IS PAN NO IN 2ND LINE SfG RANGE 1 

[401 a 

[49] -<#>./XtS 6 3 ) / C2 

[501 Z(* ]-cxt23)*(*[2]«u»-x[5]5 ) 

[511 -LOOP 
[521 G 2? 

(531 • IF RAN NO “>> l.IT MUST BE IN VD SEG 
154] a 

[551 -iapD'GZ 
[561 2(*>x[43 

[5^] -LOOP 

[50] G3 1 

[591 J«-iXl31>*<5:3)*i*-...\[T 61) 

[60 1 -l OOP 

v 


t 0 17 

::] » • GENERATES a PANPOn NurBER 

[21 •• BETuEEN .000001 AND 1 

[31 • • UITn EQuAl pROBaBIl 1 t > 

[4] J*< 71000900 '41 0000 00 

v 


RRNLIO)* 

V RNL XlAlBlDlKlKX 

t 1 1 aa COMPILES PROD UNIT INPUT PAR VALUES Gk 

[21 • • FROM GLOBAL MATRIX P& 

(3) • RIGHT ARGUMENT ■••NUMBER OF SETS 

( 4 ) a ONE SET IS REOD AS INPUT TO Fn CAlC 

(5) ■ 

[61 • D IS X BY 3 MATRIX TO INDEX fiy 

[7] a 

[0] 0-<X,3)M <X«3) 

[93 fiy* ( 27. 3"X ) ao 
[ 10 ) A- 27 3 ao 
till B~P& 

( 12 ) tfi-x-0 

(131 LOOP ’ 

t 1 4 ) a COMPILES a fik SET 
( 15 1 • 

(161 •»(27<x-x»l ) »G0 

[ 17 ] 4UI ]*PP(S SIXll 
(10 1 -LOOP 
(19] GO: 

(20! a fiy S£7S ARE LAMINATED 
[21 ) a 

122) -*(X «X1*X 1 ♦ 1 > t£NP 

(231 iVt lOCXlJ ] ]►* 

[24] x-1 

(251 -LOOP 
(26) END: 
v 

vunet 0 )* 

* Z-UMP XtAlT 

111 a a CONVERTS RIBBON GROuTh RATE FROM TO i 
( 2 1 • RIGHT ARGUMENT < • > PERCENT 

(3) a EXPLICIT RESULT ACTUAL GROUTH RATE 

[4] a 

(51 a CONVERT THICKNESS TO MM 
(61 a 

(71 r*Gkt«JlaiO 

10) a 

[91 a CALCULATE UMAX FOR THICKNESS. THEN UMAX*X 
( 10 ] a 

(111 MOO ) »X«X*0 .079li»( < l ) • 0 .5) 

c 



4. APPLICATION OF FORECASTING TECHNIQUE 

TO LARGE-AREA SILICON SHEET PROJECTION 


4.1 Projection of Silicon Sheet Technology Parameters 


The last quarterly report (3) concluded that single-ribbon 
growth and processing-technology improvement offered the 
best potential for achieving low-cost sheet-material 
objectives. Consequently., only technology-sensitive 
parameters (13 out of 27) are projected for this analysis, 
as shown in Table III. The first column contains the 


iaseline, and the remaining three our projections. 

TABLE m. Silicon Ribbon Technology Projection 

Item Near- 

No. Parameter Present Future 

Future 

Limit 

2 

Ribbon width (cm) 

2.5 

5 

10 

25 

3 

Ribbon growth rate (% of V ) 

v max' 

30 

40 

50 

60 

4 

Ribbon thickness (mm) 

0.30 

0.20 

0.15 

0.10 

5 

Yield of "cell grade" ribbon (%) 

70 

80 

90 

95 

6 

3 

Ribbon furnace cost ($ x 10 ) 

50 

25 

20 

15 

7 

Equipment life (years) 

7 

8 

10 

12 

9 

Equipment availability (%) 

70 

80 

90 

95 

14 

Number of technicians 

0.50 

0.25 

0.15 

0.10 

16 

Polysilicon cost ($Ag) 

65 

45 

30 

6 

17 

Poly yield to ribbon (%) 

80 

85 

90 

95 

18 

Services and supplies ($/day) 

30 

25 

20 

10 

20 

Energy to operate equipment (kWE) 

12 

11 

10 

3 

26 

Cell conversion efficiency (%) 

8 

10 

12 

15 
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Using the three projections as input to our production-unit 
model results in Table IV, where simulation numbers A3, 44, 
and 45 represent the "near-future," "future," and "limit" 
projections. Accordingly, we conclude at this point that 
silicon-sheet growth, as a technology, has the capability of 
providing the material at the required cost to meet a 
$500/kWE array-cost objective at a later point in time than 
represented by simulation number 44 ( $ 3 14 /kWE) , but well 
before the technology limit is reached (simulation number 
45, $42/kWE) . This projection is plotted in Fig. 8 (lower 
curve), together with another one, independently arrived at. 
As can be seen, the difference between the two is 
insignificant. 


TABLE IV. Near-Future, Future, and Limit Projection 
of Silicon Ribbon Technology 


UaULUZLQH Cl MD HC. S 11/04/75 

43 

44 

45 

RIBBONS GROUN SIMULTANEOUSLY 

1.00 

1.00 

1.00 

RIBBON UIDTH, CM 

5.00 

10.00 

25.00 

AUG YIELDED GROUTH RTE.SQ M/HR 

0.08 

0.29 

1.19 

COMBINED YIELD FACTOR 

0 .68 

0.81 

0.90 

DIRECT COST IN DOLS/SG METER 
EQUIPMENT CAPITAL RECOUERY 

7.01 

1.33 

0.22 

PERSONNEL 

36.13 

8.17 

1.79 

POLY SILICON COST 

30.84 

12.94 

1.55 

SERVICES /SUPPLIES 

14.88 

3.60 

0.61 

SUBTOTAL : 

88.86 

26.04 

j 4.16 

OVERHEAD COST IN DOLS/SQ METER 

20.16 

5.13 

1.03 

G*A EXPENSES IN DOLS'SQ METER 

10.90 : 

3.12 

0.52 

PROFIT IN DOLLARS /SO METER 

11.99 

3.43 

0.57 

TOTAL COST IN DOLS/SQ METER 

131.92 

37.72 

6.29 

DOLLARS PER KU 

1319.22 

314.36 

41.92 
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Fig. 8. Silicon ribbon technology projection. 


4.2 Association Between Technology Capability and Time 


As described in section 2.4, each technology parameter, in 
order to introduce chronology, is coupled to a point in time 
by means of an estimated probability that the projection 
will, in fact, materialize. Subsequently, the cumulative 
probability distribution function is evaluated numerically 
by computer, using a random number generator. The computer 
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program operates on a predefined 27-row by 8-column matrix. 
Although similar matrices could be developed for each year, 
only two were defined for this study (1980, 1985), shown as 
Tables V and VI. The first column of each matrix 
corresponds to the p roduc t ion- uni t parameter item number 
(see Fig. 1). The next column represents the baseline, the 
following three are the projections, and the remaining three 
the estimated probabilities. Only the latter are changed 
to simulate the progression of time. For instance, the 
second row refers to ribbon width, with values previously 
shown in Table II. Whenever probabilities are listed as 1, 
the Corresponding production-unit parameter is not subject 
to this probabilistic evaluation, 

4.3 Approach to Parameter Projection 

Application of this technique requires a technical 
understanding of what the technology is capable of 
achieving, an assessment of state-of-the-art, and, finally, 
judgment with respect to the timing of subsequent 
implementation. If silicon-sheet technology is to become 
cost-effective for -photovoltaic applications, then 
improvement of widjth capability is of prime importance. 
During the last few years, ribbon width has increased from 
1 to 2.5 cm, with 5 cm expected in the near-future. This 
factor-o £-5 Improvement is one reason for projecting 25 cm 
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TABLE V. 1980 Technology Projection and Estimated 
Probability Matrix 


Item 
No . I 

Present 

Near - 

Future 

Future 

Limit 

Probability 

i ! 

.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

2 

2.50 

5.00 

10.00 

25.00 

.05 

.65 

.30 

3 

30.00 

40.00 

50.00 

60.00 

.20 

.60 

.20 

4 

.30 

.20 

.15 

.10 

.25 

.60 

. 15 

5 

70.00 

ao.oo 

90.00 

95.00 

.40 

.50 

.10 

6 

50000.00 

25000.00 

20000.00 

15000.00 

.60 

.30 

.10 

7 

7.00 

8.00 

10.00 

12.00 

.10 

.60 

.30 

a 

.00 

10.00 

10.00 

10.00 

1.00 

1.00 

1.00 

9 

70.00 

eo.oo 

90.00 

95.00 

.40 

.50 

.10 

id 

.00 

.05 

.05 

.05 

1.00 

1.00 

1.00 

u 

.00 

25000.00 

25000.00 

25000.00 

1.00 

1.00 

1.00 

12 

.00 

.10 

.10 

.10 

1.00 

1.00 

1.00 

13 

.00 

70000.00 

20000.00 

20000.00 

1.00 

1.00 

1.00 

14 

.50 

.25 

.15 

.10 

.30 

.60 

.10 

15 

.00 

10000.00 

10000.00 

10000.00 

1.00 

1.00 

1.00 

16 • 

65.00 

45.00 

30.00 

6.00 

.60 

.30 

.10 

1 ? 

eo.oo 

85.00 

90.00 

95.00 

.35 

.60 

.05 

ie 

30.00 

25.00 

20.00 

10.00 

.15 

.50 

.35 

19 

.00 

.05 

.05 

.05 

1.00 

1.00 

1.00 

20 

12.00 

11.00 

10.00 

9.00 

.20 

.60 

.20 

21 ! 

.00 

50.00 

50.00 

50.00 

1.00 

1.00 

1.00 

22 i 

.00 

10.00 

10.00 

10.50 

1.00 

1.00 

1.00 

23 

.00 

10.00 

10.00 

10.00 

1 . 0 C 

1.00 

1.00 

24 

.00 

10.00 

10.00 

10.00 

1.00 

1.00 

1.00 

25 

.00 

168.00 

168 . GO 

168.00 

1.00 

1.00 

1.00 

26 

a. oo 

10.00 

12.00 

15.00 

.30 

.40 

.30 

27 , 

.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 


TABLE VI. 1985 Technology Projection and Estimated 
Probability Matrix 


Item 

No. 

Present 

Near - 

Future 

Future 

Limit 

I i 

Probability j 

l 

.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

2 

2.50 

5.00 

10.00 

25.00 

.02 

.28 

.70 

3 

30.00 

40.00 

50.00 

60.00 

.10 

.40 

.50 

4 

.30 

.20 

.15 

.10 

.10 

.60 

.30 

5 

70.00 

80.00 

90.00 

95.00 

.10 

.65 

.25 

6 

50000.00 

25000.00 

20000.00 

15000.00 

.05 

.70 

.25 

7 

7.00 

8.00 

10.00 

12.00 

.10 

.30 

.60 

8 

.00 

10.00 

10.00 

10.00 

1.00 

1.00 

1.00 

9 

70.00 

60.00 

90.00 

95.00 

.10 

.60 

.30 

10 

.00 

.05 

.05 

.05 

1.00 

1.00 

1.00 

11 

.00 

25000.00 

25000.00 

25000.00 

1.00 

1.00 

1.00 

12 

.00 

.10 

.10 

.10 

1.00 

1.00 

1.00 

13 

.00 

20000.00 

20000.00 

20000.00 

1.00 

1.00 

1 ,00 

14 

.50 

.25 

.15 

.10 

.10 

.60 

.30 

15 

. .00 

10000.00 

10000.00 

10000.00 

1.00 

1.00 

1.00 

16 

65.00 

45.00 

30.00 

6.00 

.10 

.50 

.40 

17 

, 80.00 

85.00 

90.00 

95.00 

.10 

t .60 

.30 

18 

30.00 

25.00 

20.00 

10.00 

.10 

.' ,41 

.50 

19 

! .00 

.05 

.05 

.05 

1.00 

1,00 

1.00 

20 

12.00 

11.00 

10.00 

9.00 

.10 

.30 

.60 

21 

.00 

50.00 

50.00 

50.00 

1.00 

1.00 

1.00 

22 

.00 

10.00 

10 .. 00 

10.00 

1.00 

1,00 

1.00 

23 

.00 

10.00 

10.00 

10.00 

1.00 

1.00 

1.00 

24 

.00 

10.00 

10.00 

10.00 

1.00 

1.00 

1.00 

25 

.00 

168.00 

168.00 

168.00 

- 1.00 

1.00 

1.00 

26 

8.00 

10.00 

12.00 

15.00 

.05 

.35 

.60 

27 

.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 
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Pack* to 

0F qS® 


as the width limit and reflecting a reasonably high degree 
of confidence (0.65 probability) that, by 1980, the width 
will be somewhere between 5 and 10 cm. Our current 
projection for 1985 is that the probability is good (0.70) 
that silicon sheets between 10 and 25 cm wide will be 
pulled . 

Doubling the growth rate from 30 to 60%, and, at the same 
time, reducing thickness to 0,10 mm, still represents a 
formidable, yet not impossible, challenge. This is also 
reflected in the estimated lower probabilities for growth 
rate- (1980: 0.2 0.6 0.2 ; 1985: 0.1 0.4 0.5), Yield (items 5 
and 17) and equipment availability (item 9) are difficult to 
project from laboratory conditions, but, from our experience 
in semiconductor manufacturing, tend to require diligent 
engineering effort, rather than technology breakthroughs , to 
improve. This is the basis for our optimistic outlook. 

Future silicon-sheet pullers are envisioned as mass-produced 
(^$20,000), special-purpose machin-es , highly reliable with 
automatic melt replenishment, and capable of near-unattended 
growth. If equipment meeting these conditions does not 
become available, perhaps because of insufficient demand for 
solar-cell products, then the outlook for photovoltaics 
should be reassessed. 
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4.4 Computer Evaluation 


With the method described in section 3, the production-unit 
computer model was iterated 600 times for each of the two 
years (1980, 1985). Each time, from a different set of 
input variables, energy- cap ac i ty cost at the level of sheet 
material was computed. Three of the 600 sets of input 
variables for 1980, which resulted in the highest 
($2492/kWE) and two lowest-cost ($203/kWE and $96/kWE) 
figures, are listed in Fig. 9. The result of calculation 
and of the subsequent statistical analysis can be seen in 
Tables VII and VIII. Similar information for 1985 is shown 
in Fig. 10 and Tables IX and X. 


tIMOM MU £tQS imf-fflflMf DATED li/M/75 


2.8 21.8 22.8 


a nf BWi aw 

1 KltSm GROHTH Mff »«/««" -Z" 2 • 1 1 2 . 28 4.91 . .OR 45 PCT 33 38 89 

4 Ritbm thickness, m 1.21 1.19 f.ti 

8 YIELD Of CELL MMLITY RIBBON, PCT - 


88 89 73 


6 HMtkfKmCg-w DOLL*** - 38321 41821 39878 

EQUIPMENT LIFE. YEARS - II. • 9.111.2 


-8 interest rate, percent - 11.1 n.i n.i 

9 EQUIPMENT MMILASIL1TY , PERCENT - 73 79 83 

personnel per shift per urchins 

11 11 NO. OF SUM - 1.18 1.18 1.18 AT f - 28111 28111 28111 

12 13 NO. OF ENGRS - 1.11 l.ll 0.16 AT t * 2HII 20110 2IIII 

14 18 NO. OF TECHN - 1.32 1.23 0.18 AT t - lllll 10IM lllll 

18 HOLY SILICON COST, DOLS/KG - 81 7 8 

17 POLY YISLO TO RIBBON , PERCENT - 81 89 84 . 

services and supplies 

18 CRUCISLE/DIE/PARTS COST PEN DAY - 21 23 17 DOLLARS 

19 POUER COST AT - 1.18 1.08 I. I S DOLLARS PER KUH 

21 ENERGY TO OPERATE EQUIPMENT - 11 10 11 KU 


21 22 Q/« - 50 81 SI PCT OF PERS* II 10 10 PCT OF RAU NATL COST 

23 GJmjS - 10 10 10 PERCENT OF DIRECT COST *OVERHEAD 


24 etOEll BEFORE TAX, PERCENT - 10 II 10 OF DC*0sH»G*A 


- 168 188 188 

28 CONVERSION EFFICIENCY, PERCENT - 11.39 12.99 13.63 

27 ENERGY DENSITY AT Am.XU/SO H PEAK - 111 


Fig. 9. 1980 list of input parameters resulting in the highest and two lowest 
energy-capacity costs. 
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TABLE VII. Three Sample (Highest and Two Lowest) 

Energy-Capacity Cost Calculations for 1980 


mmom qroun simultaneously 

l.M 

1.00 

1.00 

ribbon UlDTH, cti 

2.S2 

21.04 

22.70 

AV 6 YLD AREA GRUTN RTE,SQ N/ NR 

1.13 

1.39 

0.07 

combined yield factor 

l.M 

0.79 

0.01 

DIRECT COOT IN MLf/«0 NITER 
EQUIPMENT CAPITAL RECOVERY 

21.31 

2.47 

1.07 

PERSONNEL 

91.70 

7.97 

3.50 

POLY SILICON COST 

43.19 

3.91 

2.70 

SERVICES /SUPPLIES 

29.43 

3.12 

1 .39 

SUBTOTAL i 

109.70 

17.40 

0.01 

OVERHEAD COST IN DOLS/SO HETER 

40.00 

4.30 

1.96 

om EXPENSES IN DOLS/SO NETER 

23.40 

2.10 

1.00 

PROFIT IN DOLLARS/SO NETER 

29.80 

2.39 

1.16 

TOTAL COST IN DOLS/SO NETER 

203.01 

20.34 

13.03 

DOLLARS PER KU 

2491.97 

202.79 

95.50 


TABLE VUI. Statistical Analysis of 1980 

Energy-Capacity Cost Figures 


Parameter 

Sample size 

Maximum 

Minimum 

Range 

Mean 

V ariance 

Standard deviation 
Mean deviation 
Median 
Mode 


Value 

600 

2492 

96 

2396 

769.605 

105026 

324.078 

246.311 

725 

681 


37 41 61 


DATED 11/04/73 

i 'tiMCHs atotsN siHOLTAmoosiY - i i i 

a fflMOM USDTH, CM - 3.3 .73.6 24.1 

3 RIBBON CROUTN RATE.M/HR - 2.34 2.66 S.lf . .OR AS PCT 

4 RIBBON THICKNESS, MM 0.19 *.17 t.ll 

5 >ma Of CELL QUALITY RIBBON. ACT - 09 92 SI 

.miLr, DOLLARS - 23336 23998 23621 

7 EQUIPMENT LIFE, YEARS - 11. • 19.1 11.6 

6 INTEREST PATE. PERCENT - II.* 10.0 II. • 

9 EQUIPMENT AVAILABILITY, PERCENT 

PERSONNEL PER SHIFT PER MACHINE 
19 11 MO. OF SUPVS - I.9S 1.18 I.IS AT t 
12 13 MO. OF ENGRS 

14 18 MO. OF TECHN - 9.23 1.21 1.12 AT 9 

16 POLY SILICON COST, DOLS/KG - 84 7 

17 POtV YIELD TO RIBBON .PERCENT * 83 92 

SCTWCfS AMO SUAPLJfS 

18 CRUCIBLE /DIE /PARTS COST PER DAY - 29 22 15 DOLLARS 

19 POUER COST AT - I.9S 0.15 0.08 DOLLARS PER KUH 

21 ENERGY TO OPERATE EQUIPMENT - 19 10 10 KU 

21 22 Q/H - SB 50 59 PCT OF PERS* II 10 II PCT OF RAN NATL COST 

23 - 10 10 10 PERCENT OF DIRECT COST .OVERHEAD 

24 e*C£ir BEFORE TAX. PERCENT - 11 II 10 OF DC*0/H*G*A 

28 W iSSii - 168 168 168 

26 CONVERSION EFFICIENCY, PERCENT - 12.89 14.00 14.32 

27 ENERGY DENSITY AT AMl.KU/SQ N PEAK - 1 1 1 


81 84 88 

25110 25909 25000 
l.ll 0.11 l.ll AT B - 21010 29900 20101 

19909 10900 11000 
i 

87 


Fig. 10. 1985 list of input parameters resulting in the highest and two lowest 

energy-capacity costs. 


TABLE IX. Three Sample (Highest and Two Lowest) 

Energy-Capacity Cost Calculations for 1985 


RIBBONS CROON SIMULTANEOUSLY 

l.ll 

1.00 

1.00 

RIBBON UISTH, at 

3.29 

23.86 

24.05 

AUG YLD AREA GmSTN RTS, SO M/Wt 

0.08 

0.49 

1.86 

com I HID YIELD FACTO® 

0.73 

0.85 

1.69 

DIRECT COST IN DOLS/SQ METER 
EQUSFmNT CAPITAL RECOVERY 

7.86 

0.94 

1.49 

PERSONNEL 

SI. 81 

5*32 

2.62 

POLY SILICON COST 

32.94 

3.41 

2.29 

SERVICES /SUPPLIES 

18.29 

2.14 

l.ll 

SUBTOTAL i 

119.91 

11.81 

6.42 

OVERHEAD COST IN DOLS/SQ METER 

27.82 

2.98 

1.47 

£•* EXPENSES IN DOLS/SQ METER 

13.77 

1.46 

0.79 

PROFIT IN DOLLARS /SO METER 

15.15 

1.62 

0.67 

TOTAL COST IN DOLS/SQ METER 

166.65 

17.85 

9.54 

DOLLARS PER KU 

1292.98 

127.56 

66.66 


PRI&WAP. PAGE W 
m KX)R QUALITY 
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TABLE X. Statistical Analysis of 1985 

Energy-Capacity Cost Figures 


Parameter 

Sample size 

Maximum 

Minimum 

Range 

Mean 

Variance 

Standard deviation 
Mean deviation 
Median 
Mode 


Value 

600 

1293 

67 

1226 

413.728 

29022 

170.358 

128.373 

385 

410 


Accordingly, an energy-capaci 
projected, which is between 
median ($725/kWE). For 

$350/kWE, derived from the 
($385/kWE), and subsequent 
results, the procedure was it 
.time with a different str 
statistical analysis result 
$778/kWE and a median of $ 
were $369/kWE and $342 /kWE, r 


ty cost for 1980 of $750/ 
the mean ($770/kWE) an 
1985, the projected val 
mean ($414/kWE), the 
iterations. To confir 
erated another 600 times 
earn of random numbers 
ed, for 1980, in a me 
727/kWE. For 1985, the 
espectively, or about 10% 


kWE is 
d the 
ue is 
median 
m th e 
, each 
, Th e 
an of 
values 
love i 


than the first set, 


Figure 11 shows the frequency distributions 
sets of 600 energy-capacity costs. 1 1 is 


of these four 
noted that the 
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1985 distributions are grouped considerably tighter around 
the median, perhaps because our probability estimates have 
less variation for that later point in time. 



Energy -Capacity Cost ($/kWE) Energy-Capacity Cost <$AWE) 



Energy-Capacity Cost ($/kWE) Energy-Capacity Cost ($/kWE) 


Fig. 1 1 . Frequency distribution of calculated energy-capacity costs. 
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At JPL's suggestion, a further statistical test was per 
again using a different random number seed, primar 
confirm that a steady state had been reached. The re 
the 1985 calculation, shown in Table XI, confirme 


condition. As 

shown, 

2223 

calculations were 

pe 

(first column) 

before 

the 

APL workspace was 

f ul 

calculation was 

t e rm in 

a t e d . 

Intermediate result 

printed every 6 

0 lines. 

The 

second column shows 

t h 

the third the 

median. 

and 

the last three the 

s 

deviation, the 

mean deviatio 

n, and the remaining 

wo 


size. 


f o r me d , 
i 1 y to 
suit of 
d this 
rformed 
1 , and 
s were 
e me a n , 
tandard 
rkspace 


TABLE XI. Statistical Analysis of 1985 Energy -Capacity Cost 
Figures, Including Iterations 


Iter- 



Standard 

Mean 


alion 

Mean 

Median 

Deviation 

Deviation 

AWA 

3 

204 

269 

122 

06 

45204 

63 

369 

343 

146 

119 

44724 

123 

363 

347 

130 

102 

44244 

183 

359 

347 

128 

101 

43764 

243 

353 

335 

130 

103 

43284 

303 

355 

332 

143 

108 

42804 

363 

353 

331 

138 

105 

42324 

423 

351 

330 

141 

107 

41844 

463 

355 

334 

147 

111 

41364 

543 

360 

337 

153 

114 

40884 

603 

359 

338 

152 

114 

40404 

663 

362 

342 

155 

117 

39924 

723 

361 

340 

155 

117 

39444 

783 

361 

339 

153 

116 

38964 

843 

362 

340 

152 

116 

38484 

903 

361 

338 

155 

117 

38004 

963 

362 

340 

154 

117 

37524 

1023 

363 

339 

156 

lie 

37044 

1083 

362 

336 

156 

ne 

36564 

1143 

363 

336 

159 

120 

36084 

1203 

362 

336 

138 

120 

35604 

1263 

363 

336 

160 

121 

35124 

1323 

364 

338 

160 

121 

34644 

1303 

362 

335 

159 

121 

34164 

1443 

363 

335 

159 

121 

33484 

1503 

363 

336 

I5e 

119 

33204 

1563 

363 

336 

157 

119 

32724 

1 3 

363 

337 

157 

119 

32244 

1603 

364 

338 

157 

lit 

31 "’64 

1743 

363 

338 

156 

110 

31284 

1003 

364 

338 

156 

no 

30804 

1063 

365 

339 

156 

110 

30324 

1923 

365 

339 

157 

119 

29844 

1903 

365 

340 

156 

110 

29364 

2043 

365 

341 

156 

lie 

20084 

2103 

365 

341 

153 

lie 

28404 

2163 

365 

340 

156 

no 

27924 

2223 

365 

340 

157 

no 

27444 


62 



4.5 Analysis of Results 


An important consideration is whether or not this 
"bot toms-up" process-related analysis supports the 
projections of future energy-capacity cost arrived at by 
other means. To analyze this further, ERDA-pro j ected (4) 
"established silicon array costs" are plotted together with 
the results of this study in Fig. 12, where both a regular 
plot and a log plot are shown. 



5000 

4000 

3000 

2000 

1000 

0 


10 OK 



1 OK 


1000 


100 


Fig. 12. 






Generally speaking, the value of single-crystal silicon 
material should not exceed 30% of array cost. For the 
remainder of this decade, as can be seen from Fig. 12, our 
projection of energy-capacity cost at the silicon-sheet 
material level is somewhat less than half of the array cost 
projected by ERDA and thus within the accuracy level of this 
type of projection. From 1980 through 1985, however, this 
band narrows to a point in 1985 where projected 
energy- cap ac i ty cost at the level of sheet material 
( $ 350 /kWE) almost equals array cost ($500/kWE). 

In sum, it is concluded that compatibility and general 
agreement exist between this study and ERDA projections, 
particularly through 1980. 

5. CONCLUSIONS 

A new technology forecasting technique is being developed 
and applied to projecting the future cost- of energy at the 
level of silicon-sheet material. This technique is based 
upon the production-unit concept and deals with the 
economics of material manufacturing from a processing 
parameter, or "bottoms up," standpoint. From a baseline, or 
state-of-the-art, future technology capability is projected 
through full maturity-. 
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The concept of chronology is introduced by estimating the 
probability of meeting the objective associated with the 

• ■ ’ i ;• 

technology parameter at a stated point in time. From this 
probability density function, the cumulative probability 
distribution function is derived. The latter function is 
evaluated numerically, thus providing a set of input 
parameters to the production-unit model. Calculation, 

followed by subsequent iteration of this procedure, and 

final statistical analysis of the accumulated output from 
the model form the basis for the projected energy-capacity 
cost versus time relationship. 

Application of this technique results in the following 
outlook for large-area silicon sheets: 

o Silicon-sheet technology has the potential for 

achieving future low-cost material objectives for 
photovoltaic applications, if development 

milestones defined in this study are met. 

o 1980 and 1985 energy-capacity costs of $ 7 50/kWE 

and $350/kWE, respectively, at the level of 

silicon-sheet material, are projected. 

o This analysis confirms, from a silicon-sheet 

material standpoint, that ERDA-stated 
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energy-capacity cost objectives at the array level 
are achievable. However, there appears to be 

little, if any, margin for error. 

o Through 1980, a factor of ^2 difference exists 
between ERDA-pro j ected costs at the array level 
and our cost at the silicon-sheet material level. 
By 1985, this difference has essentially 
disappeared ($500/kWE and $350/kWE), which is an 
undesirable cost trend requiring further analysis. 
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FIFTH QUARTER ACTIVITY PLAN 


o Continue process studies on ribbon perfection. 

o Optimize 38-mm-wide ribbon growth. 

o Study influence of ribbon defects on lifet ime and 

solar-cell efficiency. 

o Continue work on comparative analysis of material- 

area throughput capability between capillary action 
shaping technique and Czochralski processing, using 
computer graphics. 

o Expand ribbon-growth computer model to address other 
material-processing steps. 
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